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_ . Abstract 

m . 

We show an auction-based algorithm to compute market equihbrium prices in a production model, 
("^ ' where consumers purchase items under separable nonlinear utility concave functions which satisfy 

^S) ' W.G.S{Weak Gross Substitutes); producers produce items with multiple linear production constraints. 

Our algorithm differs from previous approaches in that the prices are allowed to both increase and 
rS ,' decrease to handle changes in the production. This provides a tatonnement style algorithm which 

.^^ ' converges and provides a PTAS. The algorithm can also be extended to arbitrary convex production 

, regions and the Arrow-Debreu model. The convergence is dependent on the behavior of the marginal 

^^ ■ utility of the concave function. 

r '. 1 Introduction 

\^ [ The market equilibrium is a well studied problem for a general market model which includes production 

C/3 ■ constraints P'. Arrow and DebreufT had introduced a general production model for exchange markets and 

O I have shown proof of existence of equilibria. In the Arrow-Debreu model, each production schedule lies in a 

specified convex set. When the production model is constrained by positive production vectors only, convex 
programs have been obtained |55] [53] for linear cases. There has been considerable recent research on the 
>. ! complexity of computing equilibria [H [Till [TTl [HI [T51 [T71 [^ . 

OO ' To solve the market equilibrium problem, a number of approaches have been used including convex 

programming, auction-based algorithm[T2 and primal-dual[Tn] methods. 

Two techniques, the primal-dual schema and auction-based algorithms, have mainly been successful for 
models satisfying W.G.S.(weak gross substitutability). Other technique includes tatonnement processes. 
■^^ . Jain et al. [23i further generalize to a model with homothetic, quasi-concave utilities which is introduced in 

^^ ' [11] and [IS]. However, the paper is restricted in that the concave functions are assumed to be homogeneous 

^^ I of degree one. We differ in that we consider functions that are Weak Gross Substitutes. Jain et al.f5^ also 

give an explicit, polynomial sized convex program for the production planning model with linear utilities. 
These papers utilize convex programming. Codenotti et al. [5] consider gross-substitute functions with 
positive production constraints and provide approximation results via the ellipsoid method. Recent results 
that include production within the model include the work on price discrimination model [20] . Our results 
C^ , apply to convex production sets that are not constrained to be positive. 

For the production model when there is one constraint an auction-based algorithm is provided in [25^ . In 
this paper, we give an auction-based algorithm for a production model in which consumers have separable 
utilities for items that satisfy the weak gross substitutes property. Furthermore, producers have multiple 
linear production constraints. Note that we can also consider producers that sell items to maximize their 
profit as well as purchae materials or resources to produce them. In this model, note that each buyer 
chooses a subset of items that maximizes her utility and each producer chooses a feasible production plan 
that maximizes his profit at current prices. 

While auction methods have been applied to market equilibrium problems, the key aspect of many 
algorithms ([in][lS]) is that price discovery is monotone since goods are always sold out, i.e. overdemanded. 
In this model, the change in production plans during the course of the algorithm induces oversupply of 
goods. This implies that price discovery cannot be monotone. This is a critical difference from previous 
methods. Thus unlike previous auction-based algorithms for consumer and production models, we consider 
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auction algorithm which decrease price also. Since each producer has an arbitrary number of production 
constraints, updating prices affects production schedule. When producers chooses a bundle of items which 
also maximize consumer utilities, decreasing the prices is not required. However, when profitable items are 
not demanded by consumers, price decrease may be needed. Bounding the decrease in price can be done 
by small increments of the production plan along a gradient direction specifying increasing profits. The 
direction is obtained from the convex program describing the production plans. The step length is dictated 
by the behavior of the utility functions, in particular on the behavior of the marginal utility function. For 
simplicity of presentation our results are shown in the Fisher model with linear production constraints. 
They can be extended to arbitrary convex production regions and the Arrow-Debrcu model. 

This paper is organized as follow: In Section [2l we define the production market model and invariance 
of the algorithm to show the correctness and the convergence. We describe the overall idea on the auction- 
based algorithm and describe in market states and market state transitions in Section[2] Further, procedures 
are described in Section |4] and we show invariance conditions that ensure optimality in Section [5] In Section 
[SI we finally show the correctness of the algorithm in the algorithm and evaluate the complexity. In the 
Appendix, we will describe the details of the algorithm. 

2 Production Model 

We consider a market equilibrium problem with production, termed as MEP, with nonlinear utility functions 
and multiple linear production constraints. In here, we consider nonlinear utility functions satisfy W.G.S. 
property. The production model we consider has q producers, along with n consumers (traders) and m items. 
For simplicity, we assume that every item is produced in a quantity greater than or equal to e. 

Consumer i has a utility function, termed as Ui{Xi) — '^jUij(xij), and fixed initial endowment, e^. Xij 
represents the amount of allocation on item j to consumer i, Xi ^ (xii,Xi2, ...,2;^^) denotes the current 
allocation vector on the items to consumer i, and Uij{xij) : R^ — > R^ is the function representing the utility 
of item j to consumer i. Our assumption is that Uij is separable. We denote by Vij the first derivative of Ui 
w.r.t. Xij. 

We let P be the vector of prices of the items, where the j-th component pj represents the price of item 
j. Then, we can represent bang-per-buck, yi,j : ctij — Vij{xij)/pj, and let a^ be max bang-per-buck s.t. 
Ui = maxj aij . 

For producers, let Zsj represent the quantity on item j produced by producer s and sold or bought by 
producer s. Also, there are non-manufactured raw items defined as aj for all items. Producer s gains profit 
'YlijPj^aj when all items are sold out at the price P. We assume that the production schedule is constrained 

by a set of linear inequalities. Suppose producer s has Is linear constraints. Then, Vs,£ : ^ ■ Zsja^j < K^, 
where af and K^ are constants determined by the production schedules by producer s. Vs,j : Zgj G Z. 

Given a fixed vector of prices P, the following linear programs, CPi (P) and PPg (P) represent the optimal 
consumption and production schedule, respectively. 

Consumer nonlinear programming for consumer i 

m 

CP^{P): Maximize ^Ui{Xi) subject to 

Ejli XijPj < e, Vi (1) 

x,,>0 Vi,j (2) 



Production linear programming for producer s 



PPs{P): Maximize \^ pjZsj subject to 

E,^..<<^f Vs,£ (3) 

zsj>0 Vs,j (4) 



Furthermore, we assume that utihties for the items satisfy the Weak Gross Substitute property. Items 
are said to be weak gross substitutes for a buyer iff increasing the price of any item does not decrease the 
buyer's demand for other items. Similarly, items in an economy are said to be weak gross substitutes iff 
increasing the price of any item does not decrease the total demand of other items. 

Consider the consumer maximization problem CPi{P). Let Si{P) C i?™ be the set of optimal solutions 
of the program CPi{P). Consider another price vector P' > P. Items are gross substitutes for buyer i if 
and only if for all Xi E Si{P) there exists Xl G Si{P') such that pj — p' => Xij < a;^ . 

Note that since Uij is concave, Vij is a non-increasing function. The following result in [17] characterizes 
the class of separable concave gross substitutes utility functions. 

Lemma 2.1 JJ7| / Items are gross substitutes for buyer i if and only if for all j, yvij{y) is a non- decreasing 
function of the scalar y. 

e-equilibrium in the Production Model We define an e-equilibrium in the production model as follows: 

1. Sold- out condition: All produced items are sold to consumers within a factor of (1 + e). 

2. Bpb condition: The bang per buck on the items purchased by consumer i should be almost the same. 

3. Opt-prod condition: The current production plan almost maximizes the profit of each producer. 

4. Termination condition: Consumers spend all their endowment within a factor of (1 + e). 
The first three condition are referred to as conditions OPT and are detailed below: 

s is 

h-^i,] ■■ Xij > => v^j{xij) < aijPj < (1 + e)vij{xij) (6) 

/3 : Vs, j : pjZsj < (1 + f)PjZsj (7) 

where Zsj is the optimal amount of item j produced by s. Equilibrium is established together with the 
termination condition below: 

/4 : Vz : r, < ee, (8) 

3 Auction Algorithm for the Production model 

3.1 Notations and Preliminaries 

We define the demand set of consumer i as Di = {j : Vij{xij)/{{1 + e)pj) ^ cti < Vij{xij)/pj}. We work 
with a discretized price space, and at any instant a good is sold at at two level prices: pj and Pj/{1 + e). 
We let hij represent the amount of item j that consumer i buys at price pj , and ytj represent the amount of 
item j that consumer i purchases at price pj/{l + e). xtj which is the summation of hij and yij represents 
the total quantity of allocation on item j by consumer i. We also let r^ denote residual money which can 
be calculated as e^ — J2jiPj^v ~^ PjVij/i^ + ^))- 

We denote by Zg a set oi feasible production plan(schedule) of producer s i.e., Zg = {Z : A^ Z < Kg}, 
where A represents a mx i matrix whose row and column correspond to items and constraints, respectively. 
Z represents a production plan and Kg represents capacity constraints, Z,Ks G i?™. Note that Z is an 
unconstrained vector, that is, if a producer produces item j as a product, then Z may be positive; if a 
producer consumes item j as a material or resource, then Z may be negative. Let us define a profitable 
production plan of producer s a,s Zg — S3:g'nisxzez,{P"'" Z\A^ Z < Kg} at price P. 

We define a set of over-demanded items Od — {j : 3i s.t. ^^ Zgj < X^i-^ul- Moreover, we define a set 
of over-supplied items as Og — {j : ^^ Zgj > ^^ Xij}. 



3.2 Algorithm Overview 

We apply the auction mechanism to solve the problem of finding equilibrium in the linear production model. 
Note that in this application we allow the algorithm to decrease and increase prices. The market equilibrium 
in the production model is to find a price vector at which consumers maximize their utility and producers 
maximize their profit. Consumers want to purchase items at lower price, while producers want to sell items 
at higher price, creating complementary optimization constraints. In fact, the producers would wish to 
produce as much demand as the market can sustain, to maximize their profit. The algorithm we design 
allows producers to increase supply incrementally to change their production plan according to the current 
market price. 

The algorithm is iterative and its progress can be measured by the changes in its state during successive 
steps. The algorithm we design allows producers to increase supply incrementally to change their production 
plan according to the current market price. To prove convergence we will show that the producer's profit 
increases during phases of the algorithm. 

Overview of the algorithm The algorithm starts with an initialization , procedure initialize where 
the price vector is initialized, each producer is assigned an initial production schedule and each consumer 
is assigned goods, all quantities small enough. In procedure algorithm_main, the algorithm determines if 
consumers have no extra demand and producers have no way to increase their profit and consequently stops. 
The invariant conditions to be maintained at the end of every phase are conditions (O,® and ([7]). If at the 
beginning of a phase no consumer has any residual money left, the procedure can terminate satisfying all the 
conditions. Otherwise, it is determined if the consumer can outbid other consumers to acquire her desired 
items in procedure satisfy _demand and subsequently adjust the allocations of other consumers adjust_bpb 
to ensure condition (jH]). Each of these procedures may result in an increase in price. Consequently, the 
production schedule may required to be changed. This is done in procedure prod-reschedule. 

Algorithm 3.1 algorithm_main 



1: 

2: 

3: 

4: 

5: 

6: 

7: 

8: 

9: 

10: 

11 

12 

13 



initialize 

II bpb, opt-prod and sold-out conditions hold. 
while there exists consumer i with residual money do 
update bang-per-buck 

// bpb, opt-prod and sold-out conditions hold. 
satisfy-demand{i) 

II sold-out condition holds. 
adjust_bpb 

II bpb and sold-out conditions hold. 
prod_reschedule 

II bpb, opt-prod and sold-out conditions hold. 
end while 

// termination, bpb, opt-prod and sold-out conditions hold. 



Algorithm 3.2 satisfy _demand(i) 



consumer i with extra demand purchase items 
1: update a set of demanded items of consumer i 

if item j is available at lower price from consumer k then 
outbid{i, k,j,aij) 
update bang-per-buck 
else 

raise-price(j) 
end if 



7: 



The procedure, prod-reschedule, determines a feasible direction that improves the profit of a producer, 
by using a linear program. A small step of length 5 is taken along the feasible direction. The step length is 
chosen so that optimality conditions hold as shown later. The step provides a new production plan based 
on the current price vector P, and the production of items may decrease and increase as compared with 



Algorithm 3.3 adjust_bpb 



bang-pcr-buck condition holds by outbid 
1: while bang-pcr-buck not hold, i.e 3i : r^ > and 3j : onjPj < Vij{xij) do 
2: if item j is available at lower price from consumer k then 
3: outbid{i,k,j,aij) 

4: else 

5: raise_price{j) 

6; end if 
7; end while 



the previous production plan. Also, note that Zgj may be negative in that producer s consumes item j to 
produce other items; The change in the production of items leads to goods that are either oversupplied or 
over demanded. The demand and supply are then balanced in the procedures: baLod and baLos. 

Due to a change in the production plans, a good j may be reduced in quantity. Then the consumers 
holding the good h have reduced allocation and have money to, returned to them. To ensure their bang-pcr- 
buck conditions, consumers outbid other consumers in procedure baLod. After procedure baLod, procedure 
adjust-bpb is called to ensure that bang-per-buck conditions are met. Note that this may involve a rise in 
price. 

A change in production may also involve increased production of some particular good. In procedure 
baLos, the algorithm will balance over-supplied items so that there is no over-supplied item in the market. 

Procedure baLos balances the demand and supply of over-supplied items by four different procedures 
: purchase-money, transferjmoney, selLlprice and decrease-price. Procedure purchase_money is allowed 
for consumer with extra demand to acquire items. Unless consumer has residual money, procedure trans- 
fers-money will take money from other items to purchase her demand items. If the market has an item that 
is over-supplied, then producers will provide items at a lower price as shown in procedure selLlprice. Finally, 
Procedure decrease-price will be called if there is no available item at a higher price. 



Algorithm 3.4 prod_reschedule 



1: 

2: 

3: 

4 

5: 

6: 

7: 

8: 

9: 

10 

11 

12 

13 



w^hile execute lp_solver and if not at optimality, ys,j : 2;' ■ :— Zgj ± e' ^^ Zgj / q do 
if producer can increase profit then 
baLodij) : \/j e Od 
adjusLbpb 

II bpb condition holds. 
baLos : \/j e Os 

II bpb and sold-out conditions hold. 
check_profit 
else 

break 
end if 
end while 

// bpb, sold-out and opt-prod condition hold. 



Finally, the algorithm checks whether there is total profit increase in the overall process at each iteration. 
Procedure check_profit will be called to check whether consumers spend more than they did at the prior 
iteration. 

If the total profit does not increase, then changes of productions plans are voided. It is shown later that 
in this case the production plans are at near optimality. The algorithms then exits prod_reschedule. 

Choosing a value of e' We will assume that the producers will change their production plan by a factor 
of 1 + e'. We choose e' to ensure an approximation of (1 + e. 

We choose constants as follows: set ei s.t. \/i,i : Vij{xij) I {1 + e) < Vij{{l -\- ei)xij) < Vij{xij) and 
Vi, j : Vij{xij) < Vtj{{l - ei)xtj) < (1 + e)vij{xij). Let €2 = e^/e, e = J2i'^'. and let e' = min(ei,e2). 



3.3 States of the market 

The state of the market at tuxie instant t is represented as a 3-tuple {P {t) , Z (t) , X (t)) , where P{t),Z{t) 
and X(t) represent the price vector, the production plan and the current allocation (or demand) at time t, 
respectively. We will also use Pj(t), Zsj(t) and Xij{t) to denote the price of item j, the production plan for 
item j by producer s and the demand of consumer i for item j at time i, respectively. 

We consider the transitions of the state of the market at time i, {P{t), Z{t),X{t)), to the states at time 
t + 1, {P{t + 1), Z{t + 1), X{t + 1)) via an algorithmic process. 

State Transitions of the market We consider the state transitions : 

P{t + l) =fp{P{t),Z{t),X{t)) (9) 

X{t + l) =}x{P{t),Z{t),X{t)) (10) 

Z{t + l) =,fz{P{t),Z{t),X{t)) (11) 

where we design the functions fp, fx and fz to satisfy conditions ([S]) through ©on the state of the market. 
The state transition function fp,fx and fz are computed in algorithm_main. At each iteration conditions 
([5|) through ([7]) should be satisfied. Furthermore we will show that at termination, condition (|S]) is true and 
thus at termination the market is in a state of equilibrium. 

The function fp depends on a parameter P{t),Z{t) and X{t). In other words, Pj{t) will increase by 
a factor of (1 + e) unless 3i : yij{t) > 0. The procedure outbid may occur on either consumer's side or 
producer's side. Consumer i outbids consumer k to acquire item j, or producer s outbids consumer k to 
consume item j to produce item j' . 

To compute fx we invoke procedure satisfy -demand which allows consumer i with r^ > ee,; to purchase 
item j G Di at Pj{t) by acquiring item j from consumer k. If the demand and the supply does not match, 
then procedures outbid, baLos or baLod may be called so that X{t + 1) may change. 

To update the production schedule and compute fz, we invoke procedure prod_reschedule. According 
to P(i), a linear program solver returns an optimal production schedule. As mentioned before, a next 
production schedule is determined by the demand and the supply. 

3.4 Market State Transitions 

States of market in procedure satisfy .demand We now detail the change of the states of market 
as well as maintenance of invariants in procedure satisfy -demand inside procedure algorithm_main. In 
a market state, {P[t),Z{t),X{t)), procedure satisfy_demand calls either outbid or raise.price when Vj : 
^ijPj ~ '^ij{^ij)- Since there exists consumer i s.t. r^ > ee^, she will buy item k — argmaxj a^ . 

Note that before procedure satisfy -demand, Vj : a^- = Vij{xij{t))/pj{t). Then, for item j ^ fc, aijPjit) — 
Vij{xij(t)) and aijPj{t + l) = Vij{xi.j{t + l)). In the case of item k, a^kPkit) = Vikix^kit)) and a^kPhit + l) = 
(1 + e)vik{xik(t + 1)). However, aik will be updated, and finally, Vj : aijPj{t + 1) = Vij{xij{t + 1)). 

Since there is, P{t) and X{t) may be not equal P{t + 1) and X{t + 1), respectively. However, no change 
in production plan, Z{t) = Z{t + 1). 

States of market in procedure adjust-bpb We next consider changes to the states of the mar- 
ket and the invariants in procedure adjust-bpb inside procedure algorithm-main. In a market state, 
{P{t), Z{t),X(t)), procedure adjust-bpb calls either outbid or raise-price when 3j : aijPj < Vij{xij). For 
the corresponding consumer i, Vj : aijPj = Vij{xij) at the end of adjust-bpb. 

For item j s.t. aijPjit) < Vij(xij(t)), procedure outbid makes Pj{t) = pj{t + l),Xij{t) < Xij{t + 1); 
on the other hand, procedure raise-price makes (1 + e)pj(t) = pj{t + 1). The aim of this procedure is 
aijPjit + 1) = Vij{xij{t + 1)). 

Therefore, P{t) and X{t) may be not equal to P{t + 1) and X{t + 1), respectively. Since there is no 
change in production plan, Z{t) — Z{t + 1). 

States of market in procedure prod-reschedule Given {P{t), Z{t),X{t)) (the current market state) 
procedure prod-reschedule calls Ip-solver which returns optimal production plans, Z^ according to P{t) (as 
shown in line 2 in procedure prod-reschedule). Zs{t + 1) depends on Zs{t) and Zg as follows: 



'■^ ' *^^ ' ) Zsj{t) otherwise 

Since the production plan changes, allocations of consumers and production of producers might change. 
These changes may increase or decrease price of items. We describe more details in Sectional 

4 Other procedures in the algorithm 

In this section we will discuss other important and their properties. 

1. outbid(i, j, k, a) : consumer i with surplus will outbid consumer k to acquire item j. The quantity that 
is outbid of item j is determined by the utility function and the current allocation so as to maintain 
the invariancc. After we set the amount outbidden, we need to update the allocation and the surplus 
of both consumers i and k. 

2. purchase_money(i, j, iouersMpp/y) : foritemj G Oj , the procedure will check whether item j is demanded 
by consumer i with surplus. If exists, consumer i will purchase item j as much as the minimum of the 
quantity oversupplied and extra demand. Here we also need to update the allocation of consumer i 
on item j and the surplus. 

3. tT:an.sieT:jxioney{i,j,t over supply)'- in the case when there is no consumer with surplus for item j G Os, 
the procedure transfer_money will be called. That is, consumer i with demand on item j will give up 
item j' and spend the returning money to buy item j. The allocation of consumer i will be updated 
for item j and j' . 

4. selLlprice(i, J, toiiersuppZi/) '■ if there is consumer neither with surplus nor with items transferable, 
producers will offer item j at lower price. Remember that the algorithm accepts two-level price: a 
higher price and a lower price by a factor of (1 + e). Producers are forced to supply item j at a lower 
price, then either j is fulfilled or no consumer has item j at a higher price. 

5. decrease_price(j) : when procedure selLlprice fails to fulfill item j, we assure that the allocation of 
item j is at a lower price to any consumer. Then, procedure decrease_price will be called, and the 
while statement of procedure baLos will check whether procedures purchase_moeny, transfer_money 
or selLlprice might be called. 

6. raise_price(j) : after the procedure prod_reschedule, if there is a consumer with surplus then consumer 
i outbids other consumers. If consumer i is not satisfied then she will increase the price of item j by a 
factor of (1 + e). Procedure raise_price itself does not change the allocation, but affects invariance that 
should be maintained in the algorithm. As consumers acquire item j at higher price than previous, 
the producers may change their production plan. This will be checked in the program. 

5 Invar iances 

In this section we will prove the invariances that hold during the course of the algorithm. We first discuss 
the properties of two procedures. 

5.1 Procedure raise_price 

We consider the number of occurrences of procedure raise_price. 

Lemma 5.1 In procedure prod_reschedule at most one occurrence oj procedure raise_price is required. 

Proof: We claim that procedure raise_price occurs at most once per each item at each iteration in procedure 
prod_reschedule. 

Definitely, procedure baLos will not call procedure raise_price. Procedure prod_outbid allows producers 
to outbid consumers , but procedure raise_price will not be called. 



Procedure baLod will be called , and consumer i may potentially violate condition ^ because Wy {xij ) 
increases when Xij decreases. To adjust bang-per-buck condition, procedure adjust_bpb will be called. 

If Vfc : ykj ~ and a^pj < Vij{xij), then procedure raise_price will be called once. Note that atjPj = 
Vij{xij) < Vij{x[j) < (1 + t)aijPj, where x[j corresponds to the reduced amount of item j of consumer i. 
Since e' is chosen to satisfy "yij (xij ) <Vij{{\ — t')xij) < {l + e)vij{xij), atjPj <Vij{x[j) < {l + e)aijPj. After 
raise.price, i.e. p'j = {I + e)pj,aijp'j > Vij{x[j). 

Therefore, at most one occurrence of procedure raise_price will be enough in procedure prod_reschedule. D 



5.2 Procedure decrease_price 

Similarly, we consider the number of occurrences of procedure decrease_price. 

Lemma 5.2 Procedure decrease-price of each item occurs at most once per each iteration in procedure 
prod-reschedule. 

Proof: Let p and p' denote a previous price vector and a current price vector, respectively. Let z and z' 
denote a previous production plan and a current production plan, respectively. 

Suppose, for contradiction, that procedure decrease_price on item j occurs at least twice in an iteration. 
Remember that we have two additional variables ei and 62 and e' = min(ei, £2) in a paragraph in Section[3l 
It is enough to show the case when e' = 62 = e^/ J^i ^i- 

The money required to consume the oversupplied items is 

e' mi ^^j'^j' = *^^ mi "^sfPr/Y^ ^' 

s j' s j' i 

Calling decrease_price on item j returns the following amount of money to all consumers, 

i i 

For decrease_price to occur more than once, the required money must be greater than the money returned, 
i.e. 

^'^^fi'"'^^' > ((^^ + 2.) E h. + - E 2/^.>. 

Since Y.s E/ ^syPy <T.iei, 

e^ > ((e^ + 2e) E h^J + e E Vv)?! 

i i 

^ e^ > ((e + 2) E '^'^J + E V^^^P^ > Yl ^'jPj 

i i i 

=> e^ > E XrjPj > e^ 

i 

Since every item is demanded by at least one consumer, the money spent on item j , ^^ XijPj should be 
equal or at least e^, i.e. '^^XijPj > e^. This is a contradiction. D 

Ensuring sold-out condition ([5]): Procedures adjust_bpb and satisfy -demand will not violate the in- 
variance because there is no change in production plan and the amount of items sold does not decrease. 
When producer s changes a production plan from Zs{t) to Zs{t + 1), the following cases arise. 

• E.^.j(i + l)-E.2;,j(t + l)<0 

Procedure baLod resolves the first case, and the second case can be resolved in procedure baLos. That is, 
procedure baLod balances all items in Od] procedure baLos balances oversupplied items. 



1. For j G Od, if 3i : Xij > then Xij{t + 1) = Xij{t) — min(xij(i),^) shown m procedure baLod. 

2. For j G Os, if 3i : J g D^ and r^ > ee^, then consumer i wiU purchase item j by using procedure 
purchase_money (as shown in hue 2 and 3 in procedure baLos). It ensures that Pj{t) — pj{t + 1) and 
while J2s ^sj{t + l) = J2i 2;ij(i + 1), OPT(conditions sold-out® through bang-per-buck®) conditions 
are met. 

3. Otherwise, i.e. 3i s.t. j £ Dif] Os and r^ < eCi 

(a) To ensure that item j is sold out, the item is sold to consumer i s.t. j d Di. This can only happen if 
consumer i transfers money from a item that provides lower bang-per-buck, i.e consumer i spends 
money to purchase item j instead of item k s.t. Xik > /\ atk < en (in procedure transfcr_money) . 
The procedure transfer _money ensures that sold-out® condition is satisfied. 

(b) If item j is still over- supplied, then producers may offer item j at a lower-level price (remember 
there are 2 price levels). Note that the current price p^- does not change, but consumers will acquire 
item j more with the same amount of money. When ^j, Zgj (t-l- 1) — ^^ Xij (t+ 1) < e X^i ^ij (^+1)) 
item j will be sold out with the same money(in procedure selLlprice) . Then, Pj{t) = pj(t + 1) 
and J2i ^ij{t) < Si ^ijit + !)■ Conditions OPT will be satisfied. 

(c) Item j will not be sold out yet if Yls^sjit + 1) — J^t^iji^ + 1) > ^^i^ijit + 1). Then, in- 
stead of providing item j at lower-level price, producers will offer item j at a lower price than 
before i.e. pj{t -|- 1) = pj{t)/{l + e)(in procedure decrease_price). Then, Pj{t) > pj{t + 1) 
and J^i^ijit) < Y.i^ij{t + !)• If stih J^^^sjit + 1) > J^i^iji^ + !)> then conditions OPT 
may not be satisfied. However, then the iterations within the procedure baLos will repeat until 

E,^s,(t + l) = E.a:..7-(i + l)- 
Lemma 5.3 Condition (0j is satisfied at the end of procedure prod_reschedule. 

Ensuring bpb condition (l6|): 

Lemma 5.4 Condition (B^ is satisfied after procedures satisfy ^demand and adjust_bpb. 

Proof: In initialization, Vij{xij) < onjPj is true. When procedure adjust_bpb occurs in procedure algo- 
rithm_main, Vij{xij) < aijPj is always true since adjust_bpb iterates until Vj : Vij(xij) = otijPj. Also, 
(^ijPj < (1 + (-)vij{xij) is true. 

In procedure satisfy -demand, note that a^ is updated according to Vij{xij) and pj which implies 
that Vj : Vij{xij) — onjPj. For k — argmax^aij, consumer i purchases a;^^ amount of item k s.t. 
'"ik(x'^k) = ctikPk/i^ + e)- Then, v^jix^j) < aijPj is stih true. D 

Lemma 5.5 Condition (^ is satisfied at the end of procedure baLod. 

Proof: We claim that condition (jH]) holds at the end of baLod based on Lemma [STTI 

Let us consider procedure baLod where consumers balance their over-demand items. Let consumer i 

give up some amount on item j. Consumer i may potentially violate condition ([6|) because Vij{xij) increases 

when Xij decreases. 

However, as shown in Lemma |5.1[ raise_price occurs at most once at each iteration. Therefore, any price 

or allocation change does not violate condition ([6]). D 

Lemma 5.6 Condition bpb (0j is satisfied at the end of procedure baLos. 

Proof: Let producer s change his production plan, and j G Os- There are two possible situations in 
procedure baLos: l)procedure decrease_price is not called. 2)procedure decrease_price is executed. 

Procedure decrease_price is not called since consumers consume all amount of item j. Let consumer i 
have surplus, i.e. r^ > ee^, and j G Di. Since consumer i has surplus, she will buy item j as shown in 
procedure purchase _money. However, it will not violate condition ^ because consumers will only purchase 
item j within a factor of (1 -I- e) of bang-per-buck. That is, previously aijPj = Vij{xij), and xj- increases 
such that Vij(xij) < Vij{x[A < (1 + €)vij(xij). 



Although consumer i has surplus, the algorithm is allowed for consumers to purchase items by calling 
procedures transfer_money and selLlprice. Since the amount of production change is well defined, both 
procedures do not violate bpb condition. Recall that procedure transfer_money allows consumer i to buy 
item j to balance bang-per-buck. 

When all other procedures do not resolve the over-demand of item j, procedure decrease_price occurs. 
The occurrence of procedure decrease_price implies that consumers who want to buy item j do not have 
enough money. After procedure decrease_price occurs, let consumer i buy item j, and let x[j correspond to 
the new amount of item j of consumer i. 

Previously aijPj — Vij{xij) and after one occurrence of procedure decrease_price as shown in Lemma 

(1 + e)aij = Vij{xij)/p'j 
Consumer i now outbids to take item j. Then, 



Vij{xij)/{1 + e)p'^ < v,j{x'^j)/p'j < v,j{x,j)/p'^ 
(1 + e)ay/(l + e) < v,,{x'^^)/p'^ < (1 + e)a,, 
a,j < Vij{x'ij)/p'j < (1 + e)a,j 



Ensuring opt-prod condition ([7]): It happens that no producer reschedule their production plan. When 
there is no profit for producers, producers do not want to change their own production schedule. Even in 
the case, we show that producers will satisfy the following inequality as follow: 

Vs, j : PjZsj < PjZ*j < (1 + efpjZsj 

Lemma 5.7 When the procedure rolLback occurs, producer s still has 0{l + e)- approximation optimal profit 
according to the current prices. 

Proof: Let V{z,p) denote the current profit of producer s on production plan z according to the price 
vector p. Similarly, let V{z\p) be the next profit of producer s on production plan z' according to the price 
vector p. V{z,p) denotes the profit of producer s when he has the production plan z, the optimal production 
plan according to the price vector p. 

When the production plan shifts from V{z,p) to V{z',p), price may change to p' due to procedure 
decrease-price. The occurrence of procedure decrease_price may violate that producer s increases his profit. 
In procedure prod_reschedule, if producer s has no profit increase, then there is no production change. We 
will show that when procedure rolLback, producer s still guarantees his production profit is well bounded 
compared with the optimal production profit. We show this by proving that 

V{z,p)<il + 2e)Viz,p) 

Note that z, z', z are points on the poly-tope of multiple linear production constraints. By the property, 

z = z + {z' — z)/a => zp' — zp> -I- p'(z' — z)la 

where a =\z' — z\l\z — z\. 

V[z,p') = V(z,p') + {V{z\p') - V{z,p'))/a 

Let V{z\p') = {I + e')V{z,p'), 

V{z,p') = V{z,p') + e'V{z,p')/a 

Let e' /a < e, 

V{z,p')<{l + e)V{z,p) 

Note that the procedure decrease_price occurs at most once in procedure prod_reschedule which implies that 
yj:pj<{l + e)pr. 

V{z,p) <{1 + e)V{z,p') < (1 + 2e)V{z,p) 

a 
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6 Analysis of the algorithm 
6.1 Invariances 

Lemma 6.1 When the algorithm terminates conditions OPT are satisfied and the algorithm returns a 
0(1 + e)- approximate optimum. 

Proof: Condition sold-out ([5]) is true at the end of each iteration of the algorithm as shown in Lemma [531 
condition ([S]) holds at the end of prod_reschedule. 

Condition bang-per-buck ^ is satisfied at the end of algorithm as shown in Lemma I5.4[ bpb condition 
holds after calling procedures adjust.bpb and satisfy .demand. In the case of procedure prod_reschedule, as 
shown in Lemma 15.41 Lemma 15.51 and Lemma 15.61 condition ([6|) holds. 

Condition ([7]) is true at the end of prod_reschedule as shown 15.71 D 



6.2 Convergence 

Now, let us consider time complexity which guarantees that our algorithm converges in We show the con- 
vergence of each procedure before we provide time complexity. 

Remember that bidding is organized in rounds. In each round every consumer is picked once and reduces 
his surplus until r^ = 0. If there is no outbid, then procedure raise_price will occur followed by procedure 
prod_reschedule. In procedure prod_reschedule, the algorithm will balance between the demand and the 
production. 

Let iVo = logi+g , iVi = logi+g ^^^^^^ = logi+^ - and N2 = logi+^z — , where pmax = ma,Xj pj 

a.nd pmtji = miiijPj. 

Claim 6.1 After Ni * N2 rounds of bidding, either the algorithm terminates or a round robin completes. 

Proof: If any producer does not reschedule, then price rises. In the worst case, the maximum number of 
calls to raise_price is bounded by log]^_|_g , termed A^i. 

Pmin 

If producers gain profit, then at least one producer will increase her profit by a factor of (1 + e') of the 
previous profit. We can bound the number of occurrence as \og^__^^^, — , termed N2. 

After A^i * N2 rounds, procedure algorithm_main will be executed. If there is no consumer with residual 
money, then the algorithm will exit. Otherwise, the next round robin will occur. D 

Let Tob,Tis,Tf)d, and T^s denote the time taken for procedures outbid, lp_solver, baLod and baLos, 

respectively. Tis ^ 0{qm'^{m + l)L) and Tm + Tbs == mTob + nm where Tob = logi+e ( -) — ^^^^^ {vmax = 

maxy a;ij(0) and fmm — "cavnij Xij{a„iax))- Let \E\ be the number of nonzero utilities. Also, note that 
e' = min{e^/e, maxwij(0)/(l-|-e) = v{{l + e')e). 

Lemma 6.2 The time complexity of finding equilibrium in MEP is 0{Nq * Ni * N2{Tis + Tbd + Tbs))- 

Proof: Let us consider the worst case. For one round-robin, either procedure outbid or raise_price occurs. 
The number of iterations before a round robin occurs is bounded by iVi * A'^2 as shown in Claim 16.11 Note 
that between round-robins, residual money decreases by a factor of (1 + e') of total money of consumers. It 

means that total occurrence of round-robins is bounded by iVo — logj^^^, . Therefore, the total step is 

bounded by No ^ Ni * N2- 

Inside procedure prod_reschedule, four events can occur in a call of outbid. 

1. ykj becomes zero for some k. 

2. ri becomes zero. 

3. aij reduces by a factor of (1 -I- e). 
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4. Vij reaches aij in the inner while loop of algorithm main. 

The number of event ^ is bounded by the number of buyers having non-zero utilities on item j. The 
total number of ([T]) events is bounded by \E\ x A^i. The number of type ([2]) events is exactly equal to n 
in every round of bidding. At every type ([3]) event, bang-per-buck is reduced by a factor of (1 + e) which 

varies from to . In every round only one event type of (|4]) occurs for each buyer. 

Pmax Pmin 

Thus, outbid takes \E\ log]^,),^ _rnax_ ^ logj^^^ maa max _ logj^^^ ( - J — ^^^^^ as shown in [T^ . 

Pmin Pmin'^min ^6/ ^'^min 

Note that Ip.solver is called per each iterations. The function, lp_solver, takes 0{qm?{m + l)L), where 
m X I represents a matrix of production constraints and L represents an input size. Tip_soiver = 0{qm?(m + 
l)L). It takes mTob times for consumers to get money back and for the procedure decrease_procedure, and 
it takes nm times for producers to sell their produced items. Tm + Tbs = inTob + nm. 

Time complexity is 0(iVo * ^i * ^2(T,, + Tm + Tbs))- □ 

Theorem 6.1 Approximation equilibrium in the market equilibrium with production, MEP, can be deter- 
mined by a PTAS. 

Proof: At termination, termination condition ^ is true because r^ is low. Other conditions are true by 
Lemma l6. II which ensures the correctness. The time complexity result follows from Lemma 16.21 D 



7 Conclusion 

In this paper, we show an auction-based algorithm for a production model where consumers have nonlinear 
utility functions and producers have a set of linear capacity constraints. Our algorithm can also be extended 
to arbitrary convex production regions and the Arrow-Debreu model. 

References 

K. J. Arrow, H. D. Block and L. Hurwicz. On the Stability of the Competitive Equilibrium, Economet- 
rica, 27:82-109, 1959. 

K. J. Arrow and G. Debreu. Existence of an Equilibrium for a Competitive Economy. Econometrica, 
22:265-290, 1954. 

W. C. Brainard and H. E. Scarf. How to Compute Equilibrium Prices in 1891. Cowles Foundation 
Discussion Paper, 1272. 

[4] B. Codenotti and K. Vardarajan. Equilibrium for markets with Leontief Utilies. ICALP 2004 

[5] B. Codenotti, S. Pemmaraju and K. Varadarajan. A note on the Computation of Equilibria in Exchange 
Markets with Gross Substitutibility. Preprint 

[6] B. Codenotti, S. Pemmaraju and K. Varadarajan. On the Polynomial Time Computation of Equilibria 
for Certain Exchange Economics. SODA 2005 

[7] B. Codenotti, B. Mccune, S. Pemmaraju, R. Raman and K. Varadarajan. An Experimental Study of 
Different Approaches to Solve the Market Equilbrium Problem. ALENEX 2005. 

[8] B. Codenotti, B. Mccune, R. Raman, and K. Varadarajan. Does Theory Meet Practice? ESA 2005 

[9] B. Codenotti, B. Mccune, and K. Varadarajan. Market Equilibrium via the Excess Demand Function. 
STOC 2005, May 2224, 2005, Bahimore, Maryland, USA. 

[10] N. R. Devanur, C. Papadimitriou, A. Saberi, and V. Vazirani. Market Equilibrium via a Primal-Dual- 
Type Algorithm. FOCS 2002, 389-395. 



12 



[11] N. R. Devanur, and V. Vazirani. An Improve Approximation Scheme for Computing the Arrow-Debreu 
Prices for the Linear Case. FSTTCS 2003 

[12] N. R. Devanur, and V. Vazirani. The Spending Constraint Model for Market Equihbrium: Algorithmic, 
Existence and Uniqueness Results. In Proceedings of the 36th Annual ACM Symposium on the Theory 
of Computing, 2004. 

[13] B. Eaves. A Finite Algorithm for the Linear Exchange Model. Journal of Mathematical Economics, 
3:197-203, 1976. 

[14] E. Eisenberg and D. Gale. Consensus of Subjective Probabilities. The Pari-Mutuel Method, Annals of 
mathematical Statistics, 30, 165-168, 1959. 

[15] E. Eisenberg. Aggregation of utility functions. Management Sciences, 7(4):337-350, 1961. 

[16] L. Fleischer, R. Garg, S. Kapoor, R. Khandekar, and A. Saberi. A Fast and Simple Algorithm for 
Market Equilibria, WINE 2008. 

[17] R. Garg, and S. Kapoor. Auctin Algorithms for Market Equilibrium. In Proceedings of the 36th Annual 
ACM Symposium on the Theory of Computing, 2004. 

[18] R. Garg, and S. Kapoor. Auction Algorithms for Market Equilibrium: Price Roll-Backs and Path 
Auctions for the Fisher Model, WINE 2007. 

[19] R. Garg, S. Kapoor, and V. Vazirani. An Auction-Based Market Equlilibrium Algorithm for the Sepa- 
rable Gross Substitutability Case, APPROX 2001 

[20] G. Goal and V. V. Vazirani. A Perfect Price Discrimination Market Model with Production, and a 
(Rational) Convex Program for it, SACT, 2010. 

[21] K. Jain. A Polynomial Time Algorithm for Computing the Arrow-Debreu Market Equilibrium for Linear 
Utilities. Proceedings of the 45th Annual IEEE Symposium on Foundations of Computer Science, 286- 
294, 2006. 

[22] K. Jain and K. Varadarajan. "Equilibriua for Economies with Productions: Constant-Returns Tech- 
nologies and Production Planning Constraints" , SODA 2006. 

[23] K. Jain, V. Vazirani, and Y. Ye. Market equilibria for homothetic, quasi-concave utilities and economies 
of scale in production. Proceedings of the sixteenth annual ACM-SIAM symposium on Discrete algo- 
rithms, January 23-25, 2005. 

[24] K. Jain, M. Mahdian, and A. Saberi. Approximating Market Equilibrium. In Workshop on Approxima- 
tion Algorithms for Combinatorial Optimization APPROX 2003. 

[25] S. Kapoor, A. Mehta, and V. Vazirani. An Auction-Based Market Equilibrium Algorithm for a Pro- 
duction Model, Theoretical Computer Science 2007. 

[26] E. I. Nenakov, and M. E. Priniak. One algorithm for finding solutions of the Arrow-Debreu model. 
Kibernetica, 3:127-128, 1983. 

[27] D. J. Newman and M. E. Primak. Complexity of Circumsribed and Insribed Ellipsoid Methods for 
Solving Equilibrium Economics Models, Applied Mathematics and Computation, 52:223-231, 1992. 

[28] M. E. Primak. A converging algorithm for a Linear Exchange Model. Journal of mathematical Eco- 
nomics, 22:181-187,1993. 

A Appendix - Algorithm 



13 



Algorithm A.l algorithm_main 



initialize 

while there is extra demand, i.e. 3i : ri > eci do 

sat isfy -demand (i) 

adjust_bpb 

prod_reschcdulc 
end while 



Algorithm A. 2 initialize 



Vj : pj := e 

\/i : Ui :— niaxj Vij{0)/pj 

Vi : 2?i := argmaxj Vij{0)/pj 

\/s,j : Zsj :^ e/q 

\fi : Xij :— e, where j £ Vi 

ei s.t. yi,j : mina;,^ Vij{xij)/{1 + e) = Uy ((1 + £1)2;^ J 

e' :=min(ei,e2) 



Algorithm A. 3 satisfy _dcmand(z) 



Vj : ay :=v^j{x,j)/pj 

j := argmax/tti; 

if 3k : yfej > then 

outbid(i, k,j, aij/(\ + e)) 

aij := Vij[Xij)/Pj 
else 

raise_price(j) 
end if 



Algorithm A. 4 adjust.bpb 



while 3i : ri > and 3j : cnjPj < Vij{xij) do 
if 3k : ykj > then 
outbid(i, k,j,aij) 
else 

raise_price(j) 
end if 
end while 



Algorithm A. 5 outbid(i, j, fc, a) 



1 

2 
3 

4 
5 


h ■■= Vk'j 
t2 ■■= n/pj 

if {vij{aj) > apj) then 
else 


6 

7 


ts := min(5 : Vi.j{xij + 
end if 


8 
9 



i ™min(ii,t2,i3) 
hij :— hij + t 
n ■.= ri- tpj 


1 
2 


Vko ■■= Vkj - t 



S) = apj 
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Algorithm A. 6 raisc_pricc(j) 



Pi 



hij := 



Algorithm A. 7 prod_rcschcdulc 



1 


while true do 




2 
3 

4 
5 


get the optimal production plan, i.e. Vs : Zgj 
Vs,j:z',j :^Zsj±e'Y.s^s3/q 
if T,s T,jPjz',j > T,s HjPjZsj then 
baLod : Vj e Od 


— lp_solver{p) 


6 

7 


adjust_bpb 
baLos 




8 
9 


check.profit 
else 




10 


break 




11 


end if 




12 


end while 





Algorithm A. 8 baLod 



while 3j G Oj[ do 
if 3i : Xij > then 

t :== min(xy, X^i a^y - J2s ^s]) 

end if 
end while 



Algorithm A. 9 baLos 



9: 

10 
11 

12 
13 

14: 



while Vj e Os do 

if 3i : j € Vi and r^ > eci then 

purchase_money(i, j, Y,s ^sj ^ Y.i ^y ) 
else 

if Eli : j e Vi and r^ < ee^ then 

transfer_money(i, j, Y,s ^sj - J2i ^ij) 
end if 
if producers provide items at lower price, i.e. J^s ^sj ~ X^i ^ij — ^ X^i ^ij then 

selLlprice(i, j, ^^ ^sj - E» a^ij) 
else if Not enough items at lower price, i.e. Xs -^sj ~ Xi ^^ij > e' Xi '^y then 

decrease_price ( j ) 
end if 
end if 
end while 



Algorithm A. 10 purchase_money(i, j, io) 



t := mm{to,ri/pj) 

flij . llij -\- L 

n := ri - tp_ 



■j 
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Algorithm A. 11 transfcr_money(z, j, to) 



if to > and hij> > then 
t := mm{hij,p'j/pj,to) 
hij, :== hijf - tpj/p'j 

Iq . Lq - Z 

end if 



Algorithm A. 12 sell_lprice(z, j, to) 



1 


while to > do 


2 
3 
4 
5 


if to > e'hij then 

Vtj := y»j + (1 + ^')hij 

Lq . Zq 6 flij 

h^j := 


6 


else 


7 
8 
9 


Vij := Vij + to 

hij := hij -to/ {I + e'] 


10 


end if 


11 


end while 



Algorithm A. 13 decreasc_pricc(j) 



V« : /lij := (1 + e)hij + y^ 
Vi : j/ij := 

Vi : Ui :— niaxjv Viy jpy 
\/i '.Vi :— argmaxj' Viy jpj 



Algorithm A. 14 check_profit 



let z' and p' be vectors of the previous iteration. 



if Y.S Ej ^sjVi - J2j ^'sjPj < 7 = e' min, Y.j Pj^sj then 

roll_back 

break 
end if 
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